Characterizing the dynamical importance of network nodes and links 
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The largest eigenvalue of the adjacency matrix of the networks is a key quantity determining 
several important dynamical processes on complex networks. Based on this fact, we present a 
quantitative, objective characterization of the dynamical importance of network nodes and links in 
terms of their effect on the largest eigenvalue. We show how our characterization of the dynamical 
importance of nodes can be affected by degree-degree correlations and network community structure. 
We discuss how our characterization can be used to optimize techniques for controlling certain 
network dynamical processes and apply our results to real networks. 



In recent years, there has been much interest in the 
study of the structure of networks arising from real world 
systems, of dynamical processes taking place on net- 
works, and of how network structure impacts such dy- 
namics 0. Remarkably, the largest eigenvalue of the 
network adjacency matrix (which we denote A) has re- 
cently emerged as the key quantity determining many 
important properties for the study of a variety of dif- 
ferent dynamical network processes. Some examples are 
the following: (i) for a heterogeneous collection of chaotic 
and/or periodic dynamical systems coupled by a network 
of connections, the critical couphng strength |3] for the 
emergence of coherence is proportional to 1/A (ii) the 
critical disease contagion probability for the onset of an 
epidemic Q scales as 1/A; (iii) in percolation on a net- 
work, the condition for the emergence of a giant com- 
ponent also involves A In addition to these, there 
are other notable examples where A plays a similar role 

In many situations it might be desirable to control dy- 
namical processes that take place on networks. For ex- 
ample, in epidemic spreading, one would like to increase 
the threshold for epidemic transmission. In percolation, 
one might like to identify the key nodes holding the net- 
work together and protect them (e.g., in the transporta- 
tion network or the internet) or disrupt them (e.g., in the 
case of a terrorist network or pathogen protein network) . 
Such strategies would greatly benefit from a quantitative 
characterization of the effect of the removal of the dif- 
ferent nodes or edges in the network. We will define the 
dynamical importance of nodes and edges as the relative 
change in the largest eigenvalue of the network adjacency 
matrix upon their removal. This provides an objective 
quantification of the relative importance of the different 
elements of the network that could potentially be used to 
formulate control strategies for those network processes 



that are governed by the largest eigenvalue of the net- 
work adjacency matrix. We also will describe an efficient 
way to approximate the dynamical importance. 

We consider a network as a directed graph with N 
nodes, and we associate to it a iV x iV adjacency ma- 
trix whose elements Aij are positive if there is a link 
going from node i to node j with i ^ j and zero oth- 
erwise [An = 0). We denote the largest eigenvalue of 
A by A, where Au — Xu and A = Xv^ with u and v 
denoting the right and left eigenvectors of A. Accord- 
ing to Perron's theorem 0, of all the eigenvalues of A, 
the one with largest magnitude is real and positive and 
the components of the eigenvectors u and v all have the 
same sign (which we choose to be positive). It is often 
the case that A is well separated from the second largest 
eigenvalue. We define the dynamical importance of edge 
i j, lij, as the amount — AAy by which A decreases 
upon removal of the edge, normalized by A: 

. (1) 

Similarly, the dynamical importance of node k is defined 
in terms of the amount — AAfc by which A decreases upon 
removal of the node (or equivalently removal of all edges 
into and out of node k): 

AXk 

h ^ (2) 

We will now use a perturbative analysis in order to pro- 
vide approximations / to the dynamical importance / in 
terms of u and v. We first consider the importance of 
an edge lij . Let us denote the matrix before the removal 
of the edge by A and after the removal hy A + AA, the 
largest eigenvalue oi A + AA by A -I- AA and its corre- 
sponding right eigenvector by u + Au. For large matrices, 
it is reasonable to assume that the removal of a link or 
node has a small effect on the spectral properties of the 
network, so that Au and AA are small. Left multiplying 
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{A + AA){u + Au) = (A + AA)(u + Au) 



(3) 



2 



by and neglecting second order terms w-^AAAu and 
AAti-^Ait, we obtain AA — v'^ AAu/v'^u. Upon removal 



of edge i - 



^ j, the perturbation matrix is (AA);„ 
and therefore 



(4) 



We now examine the effect of removing node k. Upon its 
removal, the perturbation matrix is given by {AA)i„i — 
—Aim{Sik + Smk)- However, in this case we cannot assume 
Au is small as we did before, since Auk — —Uk (the 
left and right eigenvectors have zero fc*'* entry after the 
removal of node k). Therefore, we set Am — 6u — UkSk, 
where Cfc is the unit vector for the k component, and 
we assume 5u is small. Left multiplying Eq. © by 
and neglecting second order terms v'^AASu and AXv'^Su, 
we obtain AA = [v'^AAu — UkV^ AAck) / {v^ u — VkUk). 
Using the expression for AA, we get v'^AAu — ~2XukVk 
and Ukv'^ AAik = XukVk- Considering the network to be 
large {N 3> 1), we assume UkVk <IC v^u and obtain 

VkUk 



h 



(5) 



The commonly used eigenvector centrality of node 
k is defined as the eigenvector component Uk- Although 
closely related to it, Ik and lij take into account the 
possible asymmetry of A and are defined in such a way 
that they quantify the relative change in A upon removal 
of the node or link. For a review of other measures of 
node importance, see and references therein. 

We will now present examples of the dynamical im- 
portance of nodes in simulated and real networks. We 
consider unweighted networks (i.e., the nonzero elements 
Aij are identically one), but we emphasize that our con- 
siderations also apply to weighted networks. In consid- 
ering the simulated examples, we will try to mimic the 
often complex structure of real world networks. This 
complexity is evidenced by the fact that the degree dis- 
tribution in a large number of examples has been found 
to be highly heterogeneous (often following a power law 
in the so-called scale- free networks 0]), where the out- 
degree and in-degree are defined by df"* = ^"^^ 
= 'Yl!j=i-^ji- In the case of an undirected network 
A = A^ and d™ = = di. The 'degree distribution' 
P{d™ , is defined as the probability that a randomly 
chosen node has degrees c?*" and fi°"* [in the undirected 
case we write P{d) to denote the corresponding degree 
distribution]. Furthermore, it has been observed that the 
degrees at the ends of a randomly chosen edge often have 
positive or negative correlations (referred to as assorta- 
tive or disassortative mixing by degree pTl |. respectively). 
For example, the ratio 

p=(drrfr>e/(rf^)e, (6) 

where {■■■)e denotes an average over edges, {Qij)e = 
^'ijQ'ij/J2i,j ^ij^ is larger (smaller) than 1 in assor- 
tative (disassortative) networks. 



p = 1, = and Vi = d-", and thus Ik = 

C*4"/(Ef=i '^r*0- We wiU denote this reference 
importance by (For an undirected network this is 
equivalent to ranking nodes by their degree @.) 
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FIG. 1: Node dynamical importance Ik and I^. (solid line) 
as a function of log]^Q(ci™cifc''') for (a) an assortative network 
(open circles), and (b) a disassortative network (boxes). 

Our first examples consist of networks with positive 
and negative degree-degree correlations. We choose 
(somewhat arbitrarily) to examine networks in which 
the in- and out- degrees at each node are uncorrelated, 
P(d™,ci°"*) = Pi„(d™)Po«t(d°"*), and have the same 
distribution, Pin{d) — Pout{d) — P{d). The networks 
are generated by first prescribing a target degree se- 
quence ((i™,d°"*). In order to generate networks with 
a power law d egre e distribution, P{d) cx d"^, we will 
use, following the sequence of expected degrees 

c{i + io — l)~i/('i'~i) for the in-degrees, and a random 
permutation of this sequence for the out-degrees, where 
i — 1, . . . , iV, and c and are chosen to obtain a desired 
maximum and average degree. From these sequences, the 
adjacency matrix is constructed by setting Aij — 1 for 
i ^ j with probability proportional to d^^^d^ and zero 
otherwise [An — 0) (this is a generalization of the model 
in Chung et al. 12]). Finally, the following (based on a 
simplification of the method in Ref. is repeated un- 
til the network has the desired amount of degree-degree 
correlations as evidenced in the value of p: Two edges 
are chosen at random, say connecting node i to node j 
and node n to node m, and are replaced with two edges 
connecting node i to node m and node n to node j if 
sid^d^*^ + d^^df"* - dj,"d™* - d™d^"*) < 0, and are left 
alone otherwise. Setting s = +1 or — 1 we produce assor- 
tative or disassortative networks, respectively. 

We construct two different asymmetric networks of size 
N = 2000 following the algorithm above with 7 — 2.5, 
and c, iQ chosen such that (d) — 50 and d^ax — 350. 
For networks (i) and (ii) we used s = 4-1 and s — —1, 
respectively, until a desired value of p was reached. This 
resulted in values of p of 1.6 and 0.69 for networks (i) 
and (ii), respectively. 

In Fig. n we show on a logarithmic scale (base 10) 
our approximation to the node dynamical importance 
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Ik versus \/d^^^d^ for networks (i) (open circles) and 
(ii) (boxes), and the reference importance /° (solid line). 
We see that for the assortative networ k there is a rough 
monotonic relation between Ik and ^JcPj^^d^ , while for 
the disassortative case a functional relation ship eve n less 
clear and the nodes with low value of \f^k^^ have 
importance comparable to the highly connected nodes 
(since, due to the disassortativity, they act as bridges 
connecting the hubs). In both cases I^ and its approxi- 
mation by I^ [Eq. ISJ] are essentially the same to within 
the size of the plotted points in the figure (this will also 
apply to Fig.[2Jl. 
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FIG. 2: Node dynamical importance /fc for the two- 
community network described in the text (stars) and the un- 
correlated reference 7^ (solid line). 

Our next example is motivated by the fact that it 
is sometimes observed that real networks can be sub- 
divided into more or less well defined communities which 
have different statistics, and thus potentially different 
dynamical importance. As a simple model of such sit- 
uation, we specify a division of the nodes in the network 
into two groups of the same size, A and B, (^A\]B = 
{1,2,..., A^}, A{\B — 0), and then we construct a net- 
work following the steps above with s = -1-1, but only 
rewire the edges if all the nodes in consideration belong 
all to group A {i, j, n, m £ A). The effect of this division 
is to create a subnetwork (group A) with a correlation 
that is larger than that for the whole network. 

In Fig. |2] we show the node dynamical importance Ik 
for this network (stars) and the uncorrelated reference I^ 
(solid line). We see that the dynamical importance cap- 
tures the subdivision existing in the network. Nodes in 
the assortative region A are more dynamically important 
than nodes with the same connectivity outside 
of this region. This shows that the node dynamical im- 
portance can depend on the subdivision of the network 
into communities, and suggests that in networks with 
strong community subdivision the node dynamical im- 
portance could be weakly correlated with the degree. 

We will now consider the dynamical importance of the 
nodes in the undirected yeast protein interaction network 
[Ullil (iV = 2361), the directed Kiel University email 
network (TV = 64385), and the undirected internet 
autonomous systems (AS) network 15] {N = 21885). 
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FIG. 3: Logarithm of the dynamical importance 1^ as a func- 
tion of the logarithm of \/ d]T'd°y-^ for (a) the yeast protein 
interaction network fl^ . (b) the Kiel University email 
network [l^ . and (c) the internet (AS) network p^ . Figure 
(d) shows 7fc versus Ih for the AS network. The solid line is 
the identity. 

(Each one of these networks is an incomplete sample of 
a larger network. For the purpose of illustrating our 
ideas, we study the dynamical importance of the reported 
nodes.) The dynamical importance of the nodes in these 
three networks is shown as a function of va™(?™* in a 
double logarithmic scale (base 10) in Figs. 1 (a,b,c). The 
points were calculated from Eq. , except the rightmost 
point in Fig. 1 (b), for which our assumption VkUk <C v^u 
was not satisfied and for which we calculated Ik directly 
from the definition. Otherwise, the approximation Ik 
yielded good results, except for a relatively small bias 
towards larger values for VkUk/v^u ~ 0.1. This is illus- 
trated in Fig. 1 (d), which shows that, in spite of the 
deviation of Ik from Ik, the relationship is still mono- 
tonic, and hence does not affect the relative ranking of 
nodes. 

A striking characteristic of the three plots is that, al- 
though there is a correlation between dynamical impor- 
tance and the connectivity measured by d™d°"*, there 
are huge variations of importance among nodes of com- 
parable connectivity. In particular, for the directed 
email network [Fig.|31(b)], many poorly connected nodes 
((i™(i°"* ~ 1 — 5) have a dynamical importance vastly 
larger than some of the most connected nodes ((i'"(i°"* ~ 
10^). This suggests that, when enough information about 
the network is available, the dynamical importance of 
nodes should be a key element in the formulation of op- 
timum immunization strategies. 

We will now show how knowledge of the dynamical im- 
portance of nodes can be used to optimally reduce A in 
order to control various dynamical processes as discussed 
above. For the three networks presented above, we suc- 
cessively remove either (i) the most dynamically impor- 
tant nodes [as determined by Eq. JSJ] , (ii) the nodes with 
the highest value of d*"'^""*, or (iii) random nodes. (Af- 
ter removing a node, we recalculate the importance and 
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FIG. 4: Largest eigenvalue A(m) of the network resulting from 
removing m nodes in order of decreasing dynamical impor- 
tance (solid lines), degree (short dashed lines), and randomly 
(long dashed lines), for (a) the yeast protein interaction net- 
work [fKill, (b) the internet AS network [ij, and (c) the 
Kiel University email network |T^ . 

the degrees.) In Fig. ^ we show A(to)/A(0) as a function 
of TO, where m is the number of removed nodes and A(to) 
is the largest eigenvalue of the resulting network. We see 
that using the dynamical importance (solid hues) greatly 
improves the results over using the degree (short dashed 
lines). These two methods are, of course, much more effi- 
cient than removing nodes randomly (long dashed lines) . 

We now offer some additional general comments, (i) 
The perturbation technique used to obtain Eq. |(SJ| also 
can be used to estimate the increase of A upon addition 
of a new node and its links or the addition or change in 
the weight of one link. If we are allowed to add a new 
node with a prescribed number of in- and out-links or 
to increase the weight of a link by a specified amount, 
this approach can determine how to proceed so as to 
maximally increase A (as in a case in which one would like 



to promote synchronization '^). (ii) The perturbation 
analysis can also be applied to weighted networks that 
have some negative weight links {Aij < for some 
In particular, if the number of negative weight links is not 
too great, A is still real and positive and the components 
of u and v are still positive. In that case, Eq. shows 
that removing such a link increases A and, furthermore, 
one can use it to determine the best negative weight link 
to remove if one wishes to most increase A. (iii) Another 
use of Eq. © might be to determine from a given network 
a reduced simpler network with fewer nodes, but almost 
the same dynamics (in the sense of having almost the 
same A). This could be done by successive removal of the 
nodes of lowest Ik (as in the idea of the 'k-core' 0), and 
offers a potential tool for facilitating the understanding of 
the dynamics on a complex network, (iv) As compared 
to degree-based node ranking, the approximation Ik of 
Eq. requires computation of the eigenvectors u and v 
and hence more complete network information. 

In conclusion, we have presented an objective, quanti- 
tative measure of the dynamical importance of edges and 
nodes in a network. The dynamical importance of a node 
or edge measures how the largest eigenvalue, which con- 
trols various important dynamical processes in networks, 
changes upon removal of the given node or edge. We 
have shown how knowledge of the dynamical importance 
of nodes can be used to optimize strategies to control dy- 
namical processes dependent on the largest eigenvalue of 
the adjacency matrix of the network. 
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